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Research  interest  in  the  conversion  of  lignocellulosic  biomass  into  energy  and  fuels  through  the  pyrolysis 
process  has  increased  significantly  in  the  last  decade  as  the  necessity  for  a  renewable  source  of  carbon 
has  become  more  evident.  For  optimal  design  of  pyrolysis  reactors,  an  understanding  of  the  pyrolysis 
kinetics  of  lignocellulosic  biomass  is  of  fundamental  importance.  The  distributed  activation  energy 
model  (DAEM)  has  been  usually  used  to  describe  the  pyrolysis  kinetics  of  lignocellulosic  biomass.  In  this 
review,  we  start  with  the  derivation  of  the  DAEM.  After  an  overview  of  the  activation  energy  distribution 
and  frequency  factor  in  the  DAEM,  we  focus  on  the  numerical  calculation  and  parameter  estimation 
methods  of  the  DAEM.  Finally,  this  review  summarizes  recent  results  published  in  the  literature  for  the 
application  of  the  DAEM  to  the  pyrolysis  kinetics  of  lignocellulosic  biomass. 

©  2014  Elsevier  Ltd.  All  rights  reserved. 
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1.  Introduction 

The  pyrolysis  of  lignocellulosic  biomass  has  attracted  considerable 
attention  because  it  is  the  thermochemical  process  that  converts 
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lignocellulosic  biomass  into  liquid  fuel  (bio-oil),  which  is  easy  to  store 
and  transport  when  compared  to  raw  biomass  [1,2].  Knowledge  of  the 
kinetics  is  important  in  the  computational  fluid  dynamics  modeling  of 
biomass  pyrolysis  process  [3],  And  in  the  pyrolysis  process  of 
lignocellulosic  biomass,  the  relative  rates  of  decomposition,  cracking, 
and  repolymerization/condensation  reactions  influence  the  quantity 
and  quality  of  bio-oil  produced  as  well  as  the  long-term  storage 
stability  of  the  bio-oil  [4], 

The  distributed  activation  energy  model  (DAEM)  is  a  multiple 
reaction  model,  which  is  widely  used  in  the  pyrolysis  of 
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lignocellulosic  biomass  [5],  The  model  assumes  the  decomposition 
mechanism  take  a  large  number  of  independent,  parallel,  first- 
order  or  nth-order  reactions  with  different  activation  energies 
reflecting  variations  in  the  bond  strengths  of  species.  The  differ¬ 
ence  in  activation  energies  can  be  represented  by  a  continuous 
distribution  function  [6],  The  DAEM  is  not  only  used  to  describe 
the  pyrolysis  kinetics  of  lignocellulosic  biomass  and  its  main 
components,  but  also  to  describe  the  pyrolysis  kinetics  of  coal 
[7-9]  and  other  thermally  degradable  materials  [10-13], 

In  order  to  better  understand  the  DAEM  and  its  application  in 
the  pyrolysis  of  lignocellulosic  biomass,  an  overview  of  the 
research  and  development  in  the  DAEM  for  the  pyrolysis  of 
lignocellulosic  biomass  is  conducted.  The  derivation  of  the  DAEM 
was  first  presented.  Then,  the  activation  energy  distribution  and 
frequency  factor  in  the  DAEM  are  overviewed  and  the  numerical 
calculation  of  the  DAEM  is  discussed.  In  Section  5,  the  methods  for 
the  estimation  of  the  DAEM  kinetic  parameters  are  critically 
studied.  Finally,  the  up-to-date  research  efforts  in  the  DAEM  for 
the  pyrolysis  of  lignocellulosic  biomass  are  summarized. 


2.  Derivation  of  DAEM 


The  integration  of  Eq.  (4)  leads  to 

V  f  1  —  f1— f1  — n)/o^  exP(-Kr)dt] V<1  "V(E)d£  n#  1  ^ 

^  {  1-  /0“exp[-  f0A  exp(-£)dt]/(E>dE  n=  1 

where  a  is  the  degree  of  conversion. 

Differentiating  Eq.  (5)  with  respect  to  t,  the  expression  of 
da/dt  is 

da  exP  (-]fr)[l -O-n) foA  exP  (-ifr)df]"/n  "’/©dE  n/1 

dt  {  Jo” 71  exp  [-&-/0A  exp  (-|-)dt]/(£)d£  n  =  1 

(6) 

Most  laboratory  experiments  involving  the  pyrolysis  of  ligno¬ 
cellulosic  biomass  are  performed  at  temperatures  linearly  varying 
with  time  from  a  starting  temperature,  T0  [17],  Under  the  linear 
heating  program,  Eqs.  (5)  and  (6)  can  be  rewritten  as  Eqs.  (7)  and  (8) 

|  Wo”  [!-(! -n)frj  exp(-T,)dr]  V0  "/(E)dE  n=Al 
\  1  -  ./<T  exP  [-  Jto f  exP ( - wr) dr]/(E)d£  n  =  1 

(7) 


The  concept  of  an  activation  energy  distribution  was  originally 
proposed  by  Vand  [14]  and  later  used  by  Pitt  [15]  and  Anthony  [16], 
The  DAEM  for  lignocellulosic  biomass  pyrolysis  may  be  applied  to 
either  the  total  amount  of  volatiles  released,  or  to  the  amount  of  an 
individual  volatile  constituent.  The  description  in  this  work  follows 
the  total  amount  volatiles  released.  It  is  assumed  that  lignocellulosic 
biomass  has  many  chemical  groups  which  are  numbered  i=l,...,s, 
and  the  released  mass  fraction  for  the  ith  group  is  V,(t).  The 
contribution  of  the  decomposition  of  the  ith  group  is  assumed  to 
be  described  by  the  pseudo-nth-order  rate  equation 


d(V,/v,*)  ,, 

dt  \  Vi*  J 


=  Aj  exp 


(1) 


where  V*  is  the  total  released  mass  fraction  for  the  ith  group,  n  is 
the  reaction  order  and  k,  is  the  reaction  rate  constant,  A,  is  the 
frequency  factor,  £,■  is  the  activation  energy,  R  is  the  universal  gas 
constant,  T  is  the  temperature,  and  t  is  the  time. 

Eq.  (1)  can  be  written  in  its  integral  form 

Vi* 


j  1-  [l-O-nl/oAj  exp(-;|)df]  ‘  "  n^l 
[  l-exp[-  f„Ai  exp(-j§.)dt]  n  =  l 


The  complexity  of  lignocellulosic  biomass  pyrolysis  is  such  that 
a  continuous  distribution  f[E)  of  activation  energies  is  assumed 
where  Jjf +aE/(£)d£  describes  the  probability  that  chemical  groups 
with  a  sample  have  an  activation  energy  in  this  range  E  to  £+A£. 
The  mass  fraction  of  potential  volatile  material  with  activation 
energies  between  £  and  £+d£,  dV*,  at  a  given  time  t  is  V*J[£)d£. 

dV*  =  V*f(E)dE  (3) 


Then,  V*  and  V,  are  replaced  by  dV*  and  dV,  respectively.  Then 
Eq.  (2)  becomes 

j  V* <^1  —  [l  —  (1  —  n)  f‘A  exp(-T,)dtj1/<1  n>|/(£)d£  n*l 
|  V*{l-exp[-  ft  A  exp(-^.)dt]}/(£)d£  n=  1 

(4) 


da_|/o“?exP(-H')[1-<1”n)/r0?exP(-m-)dT]"/"  n>/(E)d£  11*  1 
dT  [Jo“£exP  [-m-  J?„^exP(-)&)dr]/©dE  n  =  1 


where  ()  is  the  heating  rate. 


3.  J\E)  and  A  in  DAEM 

A  popular  choice  for  the  activation  energy  distribution  is  the 
Gaussian  distribution  centered  at  £0  with  standard  deviation  a 


/(£)  = 


1 


OsFTn 


exp 


(£ — £0) 


2cr2 


(9) 


The  disadvantage  of  the  Gaussian  distribution  is  that  it  is 
symmetric  whereas  the  actual  reactivity  distributions  tend  to  be 
asymmetric  [18].  The  asymmetry  can  be  accounted  for  by  repla¬ 
cing  the  Gaussian  distribution  with  other  distributions,  e.g.  the 
Weibull  [19-21]  or  Logistic  [22,23]  distributions. 

The  Weibull  distribution  was  first  proposed  to  describe  chemi¬ 
cal  reactivity  distributions  for  the  pyrolysis  of  petroleum  source 
rocks  [19],  For  the  Weibull  distribution, /[£)  is  given  by 


m= 


2  (E  —  y 


exp 


E  —  y 

n 


(10) 


where  X  is  the  shape  parameter, ;/  is  the  width  parameter,  and  y  is 
the  activation  energy  threshold  (£  >  y).  The  mean  activation 
energy  £0  and  the  standard  deviation  a  of  the  Weibull  distribution 
are  defined  by 

£o  =  r+'/£^l+^)  (11) 


<i21 

where  r  is  the  Gamma  function.  Because  of  the  activation  energy 
threshold,  Eqs.  (7)  and  (8)  should  be  written  in  the  following  forms 
if  the  Weibull  distribution  is  considered 
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i-JTexp [-/lfexp(-l')dr]i(^;I)'  'exp  [-(^ 

Jd£ 

n  =  1 

(13) 

exP  (-rt)  [l -n -n)/r0  ?exP(-RT)dr]n/<1  "’KV)J  'exP  [-(V)"]d£ 

1 

/“#  exp  [-h-  si  ?  exp  ( '“p  [— pfjjjf] d£ 

n  =  1 

(14) 

The  Logistic  distribution  was  presented  in  our  previous  paper 
[22]  and  was  successfully  applied  to  describe  the  pyrolysis  kinetics 
of  cellulose  [23]  and  grape  residues  [24],  For  the  Logistic  distribu¬ 
tion,  f{E)  is  given  by 

f  E  _  *  exp  [  -  (g(E  -  ft)/V 3<p] 

\/3<7{l-|-exp  [~(jr(E—ft)/V3tJ)]j2 

where  ^  and  a  are  the  mean  value  and  standard  deviation  of  the 
Logistic  distribution. 

Capralis  et  al.  studied  the  pyrolysis  of  two  coals  and  found  that 
the  coal  pyrolysis  process  can  be  divided  into  two  steps:  the  tar 
and  light  hydrocarbon  gas  formation  during  the  primary  pyrolysis 
and  the  char  condensation,  cross-linking  reactions,  and  a  further 
gas  production  during  the  secondary  pyrolysis  [9[.  To  better 
describe  the  kinetic  behavior,  the  authors  proposed  a  double- 
Gaussian  distribution  to  present  the  reactivity  distribution  [9] 

f(E)  =  w/,  (£)  +  (!  -  w)/2(£)  (16) 

where  /i(£)  and  /2(£)  are  Gaussian  distributions,  and  w  is  the 
weight  parameter. 

There  are  three  ways  to  deal  with  the  frequency  factor  in  the 
literature.  The  conventional  method  to  use  the  DAEM  is  to  select  a 
constant  frequency  factor  value  for  all  reactions  (i.e.  1.67  x  1013  s  1 
for  coal  pyrolysis  [8],  and  101413,  1013'71  and  1013-9Os'1  for  the 
pyrolysis  of  hemicellulose,  cellulose,  and  lignin,  respectively  [25]). 
If  A  is  a  constant,  Eqs.  (7)  and  (8)  can  be  rewritten  into  the 
following  equations 

a_[1-/o”[1-(1-")?/^xp(-|=)dr]1/(1_n)/(£)d£  n*l 
{  1  -  Jo°°  exP  [  -  -p  Jr0  exp  ( -  £)  dr]/(E)d£  n  =  1 

(17) 


d«_|/o°°^exP(-i')[1-(1-n^/?„exP(-m')ciT]n/<1  "’/(£)d£  n#  1 
dT  {  Jo°>xP  [-^-fftexp  (~lrW}mdE  n  =  \ 

(18) 


In  the  two  above  equations,  the  inner  dT  integral  is 

^Texp(-i|;)dr  =  /(£,T)-/(£,To)  (19) 

where  /(£,  T)=  fg  exp  (  —  (E/RT)  )dT  is  the  temperature  integral 
[26,27], 

Second,  the  frequency  factor  in  the  DAEM  is  assumed  to  be 
connected  with  the  temperature 

A  =  A0Tm  (20) 

where  A0  is  a  constant  and  the  exponent  m  ranges  from  - 1.5  to 
2.5.  This  dependence  of  the  frequency  factor  on  the  temperature  in 
the  DAEM  was  found  to  be  used  to  describe  the  pyrolysis  kinetics 
of  lignin  [28],  peanut  shell,  pistachito  shell,  and  sunflower  shell  [29], 


In  this  case,  Eqs.  (17)  and  (18)  become 

a_|l-/o”[l-(l-<^0l'mexp(-|.)dr]1/(1"n)/(£)d£  n*  1 
{  1  -  /o°°  exP  Jt0  Tmex p  (-^)dr]/(£)d£  n  =  1 

(21) 

da _  f  JTnr^P]— ra-)[l-(l-n)rJ?oTmexP(  — CT)dT]n/<1  "V©d£  n#l 
dr  {  /o”T exP[-ifr-f TmexP(-i)dr]/(£)d£  n  =  1 

(22) 

In  Eqs.  (21)  and  (22),  the  inner  dT  integral  becomes 

£  Tmexp(--^dT  =  Im(E,T)-Im(E,T0)  (23) 

where  Im(E, T )  =  fg  Tmex p  (—(E/RT))dT  is  the  general  temperature 
integral  [30],  Here  we  especially  mention  the  temperature  integral 
or  the  general  temperature  integral  in  order  to  better  state  the 
numerical  solution  of  the  DAEM  in  the  next  section. 

Finally,  the  frequency  factor  in  the  DAEM  is  assumed  to  follow 
the  compensation  effect  (the  effect  of  an  increase  of  £0  is  partially 
or  completely  offset  by  a  ‘compensatory’  increase  in  A  [31]). 
The  relationship  between  £0  and  A  is  most  usually  represented 
[25,29,32] 

lnA  =  a+b£0  (24) 

where  £0  is  the  mean  value  of  activation  energy  distribution,  a  and 
b  are  constants.  The  value  of  a  is  usually  close  to  zero  [33]. 


4.  Numerical  calculation  of  DAEM 


Due  to  the  fact  that  there  is  an  inner  dT  integral  and  outer  d£ 
integral  in  the  DAEM,  it  is  very  difficult  to  obtain  the  exact 
analytical  solution  of  the  DAEM.  In  the  literature,  many  mathe¬ 
matical  approximate  approaches  were  pursued  to  deal  with  the 
DAEM.  Suuberg  expressed  the  temperature  integral  over  £  in  terms 
of  complementary  error  function  [34],  Du  et  al.  used  the  step 
function  approximation  technique  to  determine  the  activation 
energy  distribution  [35].  Donskoi  and  McElwain  used  high  order 
Gauss-Hermite  integration  to  evaluate  the  outer  d£  integral  [36], 
Rostami  et  al.  presented  how  it  is  possible  to  obtain  an  approx¬ 
imate  closed  form  of  the  model  by  assuming  that  as  conversion 
proceeds,  the  functional  groups  with  the  lowest  activation  ener¬ 
gies  are  lost  first  [37], 

Since  it  is  difficult  to  analytically  solve  the  DAEM,  several 
researchers  dealt  with  the  DAEM  by  means  of  numerical  techni¬ 
ques.  Braun  and  Burnham  replaced  the  upper  and  lower  limits  of 
the  outer  integral  with  £0+3<r  and  £0- 3u  and  divided  the 
activation  energy  space  into  96  equal  intervals  [38],  Gune?  and 
Gune$  used  5001<Jmol_1  as  the  upper  limit  of  the  outer  d£ 
integral  and  50  integral  intervals  [39], 

For  purposes  of  numerical  calculation  of  the  DAEM,  either 
general  purposed  mathematical  software  or  a  computer  program 
developed  in  any  programming  language  is  used.  In  this  work,  we 
developed  a  program  written  in  the  Mathematica  software  system 
to  numerically  solve  the  DAEM. 

Either  the  temperature  integral  or  the  general  temperature 
integral  has  no  exact  analytical  solution  [40],  but  they  can  be 
numerically  solved.  In  the  Mathematica  software  system,  the 
temperature  integral  and  general  temperature  integral  can  be 
expressed  in  the  following  forms 


f(E,T)  =  | 


Gamma 


£  ’ 
RT 


(25) 
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7E\m  +  1 

£1 

-  Gamma 

—  1  —  m,— 

W 

rt\ 

For  the  outer  d£  integral,  we  use  the  ‘NIntegrate’  function  of  the 
Mathematica  software  system  with  the  integral  limits  0  and 
£o+30cr. 

To  illustrate  the  proceeding  method,  the  following  theoretically 
simulated  DAEM  processes  with  Gaussian  distribution  are  consid¬ 
ered.  The  kinetic  parameter  values  of  the  first  process  are  as 


Fig.  1.  The  simulated  a-T  and  da/dT-T  curves  of  the  Gaussian  DAEM  processes  with 
the  following  parameters:  A— 1.67  x  1013  s  ;,  E0— 228 1<J  mol  1 , /2 —  10 1<  min-1,  n=1.0, 
and  a =0.5%  E0,  2.5%  E0,  5.0%  E0,  7.5%  Eg  and  10.0%  E0. 


Fig.  2.  The  simulated  a-T  and  da/dT-T  curves  of  the  Gaussian  DAEM  processes 
with  the  following  parameters:  A=1.67  x  1013  s  Eg— 223  kj  mol-1,  [>— 10  Kmin-1, 
(7=5.0%  E0,  and  n=1.0, 1.5,  2.0,  2.5  and  3.0. 


follows:  A=1.67  x  1013  s-1,  £o=228  kj  mol-1,  /3=10  Kmin-1,  n=l 
and  <7=0.5 %£0,  2.5%  £0,  5.0%  £0,  7.5%  £0  and  10.0%  £0.  The  values  of 
A,  E0  and  /?  were  taken  from  the  literature  [8],  which  were  used  to 
describe  the  pyrolysis  kinetics  of  coal.  The  parameter  values  in  the 
second  process  are  identical  to  the  first  example  except  for 
<7=5.0%  £0  and  n  =  1.0,  1.5,  2.0,  2.5,  and  3.0  (those  n  values  were 
used  to  describe  the  pyrolysis  kinetics  of  biomass  in  the  literature 
[41]).  Figs.  1  and  2  show  the  a—T  and  da/dT-T  curves  for  the  first 
and  second  processes.  From  Figs.  1  and  2,  it  can  be  seen  that 
increasing  the  values  of  <7  makes  the  reaction  profile  broader  and 
more  symmetric  and  that  increasing  the  reaction  order  makes 
the  reaction  profile  broader  and  makes  it  more  skewed  to  high 
temperature. 

5.  Parameter  estimation  methods 

The  DAEM  has  some  mathematical  difficulties  in  parameter 
estimation  because  of  the  complex  structure  of  the  DAEM  equation 
[28].  The  parameter  estimation  methods  for  the  DAEM  can  be 
classified  as  distribution-free  and  distribution-fitting  methods. 
Miura  and  his  colleague  [42,43]  proposed  the  differential  and 
integral  methods  to  estimate  f(E)  in  the  DAEM  without  previous 
assumptions  for  f[E )  in  1995  and  1998,  respectively.  In  further 
discussion,  they  are  referred  to  as  the  Miura  differential  method 
and  the  Miura-Maki  integral  method  [6],  which  are  distribution- 
free  methods.  The  previous  assumption  for  f{E)  is  required  in  the 
use  of  the  distribution-fitting  methods,  which  are  based  on  some 
direct  search  optimization  methods. 

5.3.  Distribution-free  methods 

The  assumptions  and  detailed  derivation  processes  of  the 
Miura  differential  method  and  the  Miura-Maki  integral  method 
can  be  found  in  the  literature  [42,43],  The  procedures  for  calculat¬ 
ing  f[E )  by  using  these  distribution-free  methods  are  compared  in 
Table  1.  From  Table  1,  it  can  be  obtained  that  the  a—T  data  at 
different  heating  rates  are  required  in  the  use  of  the  Miura-Maki 
integral  method,  while  the  Miura  differential  method  employs 
instantaneous  a-T  and  da/dT-T  data.  Therefore,  for  the  Miura 
differential  method,  it  is  more  sensitive  to  experimental  noise. 

To  evaluate  the  validity  of  the  Miura  differential  method  and 
the  Miura-Maki  integral  method,  it  is  examine  to  see  if  f[E ) 
determined  from  the  a-T  and  da/dT-T  data  constructed  from 
known  sets  of  DAEM  parameters  is  consistent  with  the  assumed 


Table  2 

DAEM  parameters  of  cellulose,  hemicellulose  and  lignin  presented  in  the  literature  [25]. 


Material 

£0  (kj  mol  *) 

o  (kj  mol  ') 

A(s-1) 

Hemicellulose 

175.6 

4.3 

io1452 

Cellulose 

185.4 

3.5 

io13-64 

Lignin 

195.4 

32.0 

io13-98 

Table  1 

Comparison  between  the  Miura  differential  method  and  the  Miura-Maki  integral  method  to  estimate  /[£). 


Step  The  Miura  differential  method 


The  Miura-Maki  integral  method 


1  Obtain  a  versus  7  and  da/dT  versus  7  relationships  at  three  heating  rates  at  least  Obtain  a  versus  7  relationships  at  three  heating  rates  at  least 

2  Calculate  the  values  of  ln[(da/d7)/(l  -  a)]  at  selected  a  values  for  different  heating  Calculate  the  values  of  ln(/?/71 2 3 4)  at  selected  a  values  for  different  heating 

rates  rates 

3  Plot  ln[(da/d7)/(l  -  a)]  versus  1/7  at  selected  a  values,  and  determine  £  Plot  ln(/?/72)  versus  1/7  at  selecteda  values,  and  determine  £ 

4  Plot  a  versus  £,  and  differentiate  a  versus  £  relationship  to  obtain  f[E)  Plot  a  versus  £,  and  differentiate  a  versus  £  relationship  to  obtain  f [£) 


JOUI  f>[  /  3  |0U1  /  (3)g 
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a  a 


a 

Fig.  3.  The  activation  energy  values  as  a  function  of  conversion  for  the  pyrolysis  of  hemicellulose,  cellulose  and  lignin. 


El  kJ  mol'1 


Fig.  4.  Comparison  of  the  assumed  f(E)  and  the  j{E )  estimated  by  the  Miura  differential  method  and  the  Miura-Maki  integral  method. 
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Table  3 

Different  sets  of  kinetic  parameters  of  the  first-order  Gaussian  DAEM  for  cellulose 
pyrolysis. 


Set  of  kinetic 
parameters 

A  (s  ’) 

E0 

(kj  mol-1) 

G 

(kj  mol-1) 

SSR 

1 

5.5472  x  1011 

161.8438 

4.3217 

0.012030 

2 

6.3165  x  1013 

185.1297 

4.5979 

0.002649 

3 

1.0107  x  1014 

187.4462 

4.5818 

0.002569 

4 

9.9504  x  1014 

198.7059 

5.2204 

0.003468 

5 

1.0004  x  1016 

210.0828 

5.8349 

0.004290 

6 

9.8514  x  1017 

232.7228 

6.9938 

0.005719 

one.  For  this  purpose,  some  simulated  DAEM  processes  are 
considered. 

In  the  literature  [25],  the  first-order  Gaussian  DAEM  was  used 
to  successfully  describe  the  pyrolysis  kinetics  of  hemicellulose, 
cellulose  and  lignin,  which  are  three  main  components  of  ligno- 
cellulosic  biomass.  Their  parameters  are  employed  to  construct  the 
simulated  DAEM  processes,  as  shown  in  Table  2. 

Fig.  3  shows  the  activation  energy  values  as  a  function  of 
conversion  for  the  pyrolysis  of  hemicellulose,  cellulose  and  lignin. 
From  Fig.  3,  it  can  be  observed  that  the  activation  energies 
significantly  increase  with  increasing  conversion.  Differentiating 


Fig.  5.  Multiplicity  in  kinetic  parameters  of  the  first-order  Gaussian  DAEM  for  cellulose  pyrolysis. 
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the  a  versus  E  relationships  gives  fiE)  curves,  which  are  shown  in 
Fig.  4  with  the  pre-assumed  Gaussian  distributions.  The  estimated 
activation  energy  distributions  are  close  with  the  assumed  dis¬ 
tribution  for  lignin.  However,  the  errors  involved  in  the  activation 
energy  distributions  obtained  from  the  Miura  differential  method 
and  the  Miura-Maki  integral  method  are  significant  for  hemicel- 
lulose  and  cellulose.  Therefore,  the  validity  of  the  Miura  differ¬ 
ential  method  and  the  Miura-Maki  integral  method  for  the 
estimation  of  the  DAEM  kinetic  parameters  should  be  system¬ 
atically  investigated. 

5.2.  Distribution-fitting  methods 

The  conventional  procedures  of  the  distribution-fitting  meth¬ 
ods  for  the  determination  of  the  kinetic  parameters  of  DAEM 
include  four  steps 

(1)  A  is  assumed  to  be  an  unknown  constant,  or  is  directly 
assigned  to  a  given  constant,  or  A  is  equaled  to  AoT”1; 

(2)  f[E)  is  assumed  to  be  a  Gaussian  or  another  distribution; 

(3)  The  following  objective  function  is  constructed 

SSR  =  («exp,i  —  Oxal.i)2  (27) 

i  =  1 

where  aexP  and  acal  are  experimental  and  calculated  values  of 
conversion,  respectively; 

(4)  A  and  fiE)  are  obtained  by  minimizing  the  objective  function  (27). 

In  the  literature,  the  grid  search  method  [19,44],  genetic  algo¬ 
rithm  [45],  pattern  search  method  [46],  Multistart  algorithm  [47], 
Hook-Jeeves  method  [48],  and  simulated  annealing  optimization 
method  [28]  were  usually  used  for  minimizing  Eq.  (27). 

The  disadvantage  of  the  distribution-fitting  methods  is  that 
different  sets  of  DAEM  kinetic  parameters  provide  an  approxi¬ 
mately  equally  good  fit  to  experimental  data  at  a  single  heating 
rate  [25,49],  To  demonstrate  the  existence  of  multiplicity  in  DAEM 
kinetic  parameters,  the  experimental  data  (52  data  points)  of 
cellulose  pyrolysis  at  a  heating  rate  of  10  I<  min-1  is  used.  The 
Gaussian  distribution  of  first-order  reactions  is  employed  to 
analyze  the  experimental  data  of  cellulose  pyrolysis.  Parameter 
estimation  is  performed  by  means  of  the  pattern  search  method. 
Several  sets  of  parameters  is  obtained  with  different  initial  guesses 
and  listed  in  Table  3.  The  results  from  the  first-order  Gaussian 
DAEMs  with  different  sets  of  kinetic  parameters  are  shown  in 
Fig.  5,  where  it  can  be  seen  that  the  DAEMs  compare  very  well 
with  the  experimental  data.  Fig.  6  shows  the  relationship  between 
the  multiple  sets  of  the  mean  value  ( E0 )  of  the  activation  energy 
distribution  and  the  corresponding  frequency  factor  (A).  A  perfect 
linear  relationship  can  be  obtained  between  InA  and  E0 

lnA  =  28.61335 +  4.92587E0  (28) 

where  A  and  E0  are  expressed  in  s-1  and  kj  mol-1,  respectively. 
This  kind  multiplicity  of  the  DAEM  kinetic  parameters  was  studied 
in  detail  in  the  literature  [12,49].  Therefore,  the  use  of  the  method 
involving  nonlinear  least  squares  optimization  of  multiple  heating 
rate  data  is  expected  to  result  in  unique  kinetic  parameters  [50], 


6.  DAEM  for  lignocellulosic  biomass  pyrolysis  kinetics 

Lignocellulosic  biomass  consists  of  three  major  components 
(i.e.  cellulose,  hemicellulose,  and  lignin)  and  also  contains  a  variety 
of  minor  components  (i.e.  inorganic  matters)  [51],  Cellulose, 
hemicellulose,  and  lignin  are  all  long-chain  biopolymers  [52,53], 
There  are  a  lot  of  reactions  involved  in  the  thermal  decomposition 
of  cellulose,  hemicellulose,  and  lignin.  Therefore,  it  is  not  accurate 


Fig.  6.  The  relationship  between  E0  and  InA  for  multiple  sets  of  kinetic  parameters. 
The  solid  line  represents  the  linear  fitting.  (A  is  expressed  in  s  1 ). 


enough  to  describe  the  entire  pyrolysis  processes  of  them  with  a 
single  reaction  [54],  In  fact,  the  pyrolysis  of  a  biomass  component 
requires  a  range  of  activation  energies  that  vary  with  the  size  of  the 
polymer  and  the  location  of  the  bonds  that  are  being  broken  [55], 
Therefore,  the  DAEM  has  been  used  previously  to  describe  the 
pyrolysis  of  biomass  components  [23,55-60], 

In  general,  the  conversion  rate  curve  (i.e.  the  da/dT-  T  curve)  of  the 
pyrolysis  of  lignocellulosic  biomass  shows  a  peak  and  a  peak  shoulder 
on  the  left  side  of  the  peak  [5[.  However,  according  to  the  numerical 
results  of  the  DAEM  processes  [29,39,41,61,62],  the  conversion  rate 
curves  of  the  DAEM  with  a  certain  activation  energy  distribution  (i.e. 
Gaussian,  Weibull,  Logistic  distributions)  and  certain  frequency  factor 
(i.e.  A= constant,  A=A0Tn,  A  =  ea+bE°)  show  one  single  peak  and  no 
peak  shoulder.  Thus,  it  is  impossible  to  describe  the  entire  pyrolysis 
process  of  lignocellulosic  biomass  with  such  the  DAEM.  Some 
researchers  reported  that  they  used  the  DAEM  to  describe  the 
pyrolysis  kinetics  of  lignocellulosic  biomass  based  on  the  Miura-Maki 
integral  method  [63-67],  However,  they  only  obtained  some  values  of 
Ea  and  A„  or  the  ranges  of  E  and  A  in  a  certain  range  of  conversion.  j[E) 
was  supposed  to  be  given  by  numerical  differentiating  the  a  versus  Ea 
relationship.  However,  sometimes  there  were  several  points  in  their 
obtained  a—Ea  plot,  sometimes  a  varied  with  Ea  dramatically.  In  those 
situations,  it  is  difficult  to  get  the  derivatives  of  the  a  versus  Ea 
relationship. 

Since  the  pyrolysis  kinetics  of  biomass  components  can  be 
described  well  by  using  the  single  DAEM,  some  comprehensive 
models  consisting  of  two  or  three  parallel  DAEM  reactions  were 
proposed  to  describe  the  overall  decomposition  kinetics  of  lignocellu¬ 
losic  biomass  [25,54,55,68],  In  those  models,  the  lignocellulosic 
biomass  was  regarded  as  a  sum  of  several  pseudocomponents.  It 
was  further  assumed  that  there  was  no  interaction  among  these 
pseudocomponents  and  the  pyrolysis  of  each  pseudocomponent 
would  be  described  by  using  a  single  first-order  Gaussian  DAEM. 


a(t)  = 


IfjjC, ■«,■(!-) 

IftiG 


(29) 


d«_Zfi1ci(dai/dt) 
df~  1“  ,G 


ai(t)  =1  — 


ysrexp[-/„'Aexp(- wr) 


dt- 


( E-Eoif 
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d  E 


(31) 


J.  Cai  et  al.  /  Renewable  and  Sustainable  Energy  Reviews  36  (2014)  236-246 


243 


Table  4 

Summary  of  DAEM  for  lignocellulosic  biomass  pyrolysis  kinetics. 


Biomass  Experimental  method  and  condition  Model 


Parameter  Kinetic  parameters  References 

estimation  method 


Jerusalem 
artichoke  tubers 

Thermogravimetric  analysis— linear  heating: 

5, 10,  20  and  30  K  min-1;  inert  atmosphere: 
nitrogen 

The  first-order  DAEM 

The  Miura-Maki 
integral  method 

E:  152-232.5  kj  mor1 
A:  2.93x10"- 
7.97  x  1023  s_1 

Pine  pellets 

Thermogravimetric  analysis— linear  heating: 
10, 15,  and  20  K  min-1;  inert  atmosphere: 
nitrogen 

The  first-order  DAEM 

The  Miura-Maki 
integral  method 

E:  160-270  kjmol-1 

A:  101,-1016s_1 

Sugarcane  trash 
sugarcane 
bagasse  pine 
sawdust 

Cedar  sawdust 
coffee  residue 

Thermogravimetric  analysis  (sugarcane  trash, 
sugarcane  bagasse,  pine  sawdust,  cedar 
sawdust,  coffee  residue)— linear  heating: 
0.0833-1  Ks-1;  inert  atmosphere:  nitrogen 

The  first-order  DAEM 

The  Miura-Maki 
integral  method 

Sugarcane  trash 

E:  110-320  kj  mor1 
Sugarcane  bagasse 

E:  110-320  kj  mor1 
Pine  sawdust 

E:  200-350  kj  mor1 
Cedar  sawdust 

E:  120-315  kj  mor1 

Xylan 

Thermogravimetric  analysis— linear  heating: 

The  first-order  Gaussian 

Distribution-fitting 

Xylan— 

Microcrystalline 

cellulose 

5  Kmin-1;  inert  atmosphere:  helium 

DAEM 

method:  pattern 
search  method 

£0— 178.311  kj  mor1, 
rr  =  5.848  kj  mor1, 

4  =  1012-702  s-1 
Cellulose— 

£0  =  210.037kjmor1, 
rr=0.944  kj  mor1, 

4  =  1013-904s-1 

Indulin  AT  lignin 
Lignoboost™ 
lignin 

Acetocell  lignin 

Thermogravimetric  analysis— linear  heating: 

5, 10,  and  15  K  min-1  for  indulin  AT  lignin,  5, 
15,  and  40  K  min-1  for  Lignoboost™  and 
acetocell  lignin;  inert  atmosphere:  nitrogen 

The  first-order  Gaussian  DAEM 

The  Miura-Maki 
integral  method 

Indulin  AT  lignin: 
£o=192  kj  mor1, 
o=43  kj  mor1, 
A=1013s-1 
Lignoboost™  lignin: 
£0=193  kj  mor1, 
o=37  kj  mor1, 
A=10,3s_1 

Acetocell  lignin: 
£0=193  kj  mor1, 
o=47  kj  mor1, 
A=1013s-1 

Alcell  lignin 

Kraft  lignin 

Thermogravimetric  analysis— linear  heating: 

5, 10  and  15  K  min-1;  inert  atmosphere: 
nitrogen 

Fixed-bed  reactor— linear  heating  rates:  5, 10 
and  15  Kmin-1;  inert  atmosphere:  helium 

The  first-order  DAEM 

Miura-Maki 
integral  method 

Alcell  lignin— £:  129- 
361  kj  mor1.  A: 

6.2  x  10"-9.2  x  1022  s 
£o  =  254.3  kj  mor1, 
o=29.0  kj  mor1 

Kraft  lignin— £:  80- 
158  kj  mor1.  A: 

3.3  x  107-1.84  x  109  s~3 
£0= 123.0  kj  mor 
o=11.0  kj  mor1 

Alcell  lignin 

Thermogravimetric  analysis— linear  heating: 

5, 10,  and  15  I<  min-1;  inert  atmosphere: 
nitrogen 

The  first-order  Gaussian  DAEM  with 
A=A0Tm 

Distribution-fitting 
method:  simulated 
annealing 
optimization 
method 

£0  =  135kjmor1 
o=33  1<J  mor1 
m= 2.8 

A0  =  92670min-1  K  1 

Hemicellulose 

Cellulose 

Lignin 

Wood 

Thermogravimetric  analysis— linear  heating: 

5  Kmin-1;  inert  atmosphere:  nitrogen 

The  first-order  Gaussian  DAEM 

Distribution-fitting 

method 

Hemicellulose— 

£o  =  132.9  kj  mor1, 
o=2.4  kj  mor1, 

A  =  1.25x  10los-' 

Cellulose— 

£0= 175.6  kj  mol-1, 
c 7=0.4  kj  mol-1, 

A  =  3.41  x  1012s-1 
Lignin— 

£0= 101.0  kj  mol-1, 
<t=11.3  kj  mol-1, 

A  =  2.22  x  106  s-1 
Wood— 

£0= 123.52  kj  mol-1, 
<7=5.97  kj  mol-1, 

A  =  1.57  x  108  s-1 

Fresh  water  algae  Thermogravimetric  analysis— linear  heating:  The  nth-order  Gaussian  DAEM  Distribution-fitting  5Kmin-1: 

5  and  10  K  min-1;  inter  atmosphere  method:  Multistart  £0=189.15  kj  mol-1, 

algorithm  <7=14.73  kj  mol-1, 

A  =  1.1291  x  1016s-1, 
n= 7.88 
10  K  min-1: 

£0= 190.02  kj  mol-1, 
<7=14.73  kj  mol-1, 


[66] 

[67] 

[63] 


[55] 


[58] 


[57] 


[28] 


[56] 


[47] 
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Table  4  (continued ) 


Biomass  Experimental  method  and  condition  Model  Parameter  Kinetic  parameters  References 

estimation  method 


4=2.6021  x  1016s-\ 
n=6.63 


Cellulose 

Thermogravimetric  analysis— linear  heating: 

5,  25,  and  50  K  min-1;  inert  atmosphere: 
nitrogen 

The  nth-order  Logistic  DAEM 

Distribution-fitting 
method:  pattern 
search  method 

^  =  258.5718  kjmor1 

(7=2.6601  kjmor1 

4  =  1.6218  xl017s-' 
n  =  1.1101 

[23] 

Corn  stalk  skin 

Thermogravimetric  analysis— Linear  heating: 
20,  35,  50,  75,  and  100  Kmin-1;  Inert 
atmosphere:  nitrogen 

The  first-order  DAEM 

The  Miura-Maki 
integral  method 

E:  62-169 kjmor1 

4:  10,o-8-102s-5s-' 

1651 

Rice  straw 

Rice  husk 

Corncob 

Thermogravimetric  analysis— linear  heating: 

2,  5  and  10  I<  min-1;  inert  atmosphere: 
nitrogen 

The  first-order  DAEM 

The  Miura-Maki 
integral  method 

Rice  straw— E:  118— 

208  kj  mol-1 
/( E)  peaks  at 

170  kj  mol-1 

Rice  husk— E:  ISO- 
224  kj  mol-1 
/(E)  peaks  at 

174  kj  mol-1 

Corncob— E:  MO- 
229  kj  mol-1 
/(E)  peaks  at 

183  kj  mol-1 

[64] 

Sawdust 

Thermogravimetric  analysis— linear  heating: 

5, 10, 15,  and  20  K  min-1;  atmosphere: 
hydrogen  or  syngas  (H2:  68.1%,  CO:  30.1%,  Q- 
C4:  0.9%,  C02:  0.9%) 

The  first-order  DAEM 

The  Miura-Maki 
integral  method 

E:  161.9-202.3  kjmor1 
4:  4.8  x  105- 
4.5xl013s-> 

[69] 

Wheat  straw 
Winter  barley 
straw 

Oat  straw 

Thermogravimetric  analysis— linear  heating: 

11,  22,  and  47  I<  min-1;  inert  atmosphere: 
nitrogen 

Two  parallel  first-order  Gaussian 
DAEM  reactions 

-dm/dt  =  ci  dai/dt+c2da2/dt 
where  ax  and  a2  are  described  by 
the  first-order  Gaussian  DAEM 

Distribution-fitting 

method 

See  Table  4.1 

1681 

Corn  stalk 

Rice  husk 

Sorghum  straw 
Wheat  straw 

Thermogravimetric  analysis— linear  heating: 

4  and  40  K  min-1  (corn  stalk,  rice  husk, 
sorghum  straw,  and  wheat  straw);  stepwise 
program:  20  K  min-1  with  steps  at  250,  300, 
350,  and  400  °C,  20  K  min-1  with  steps  at  225, 
275,  325,  and  375  °C  for  corn  stalk,  sorghum 
straw  and  wheat  straw;  20  K  min-1  with 
steps  at  225,  275,  325,  375,  and  425  °C, 

20  K  min-1  with  steps  at  200,  250,  300,  350, 
and  400  °C  for  rice  husk;  Inert  atmosphere: 
argon 

Three  parallel  first-order  DAEM 
reactions 

-  dm/dt  =  Y,i  =  i  cidai/dt 
where  Cj  is  the  amount  of  volatiles 
formed  from  a  unit  mass  of 
pseudocomponent  j ;  aj  is  described 
by  the  first-order  Gaussian  DAEM. 

Distribution-fitting 

method 

See  Table  4.2 

[251 

Corn  stover 

Cotton  stalk 

Palm  oil  husk 

Pine  wood 

Red  oak 

Sugarcane  bagasse 
Switchgrass 

Wheat  straw 

Thermogravimetric  analysis— linear  heating: 

5  and  20  K  min-1;  inert  atmosphere:  helium 

Three  parallel  first-order  DAEM 
reactions 

-dm/dt  =  2?=iQdaj/dt 

where  Cj  is  the  amount  of  volatiles 
formed  from  a  unit  mass  of 
pseudocomponent  j ;  aj  is  described 
by  the  first-order  Gaussian  DAEM 

Distribution-fitting 
method:  pattern 
search  method 

See  Table  4.3 

[55] 

Table  4.1 

Kinetic  parameters  of  lignocellulosic  biomass  pyrolysis  in  the  literature  [68], 


Wheat  straw 

Winter  barley  straw 

Oat  straw 

£oi  (kj  mol  ]) 

167.3 

167.3 

167.3 

Eo2  (kj  mor1) 

225.7 

225.7 

225.7 

o'!  (kj  mol-1) 

1.8 

1.8 

1.8 

o2  (kj  mol-1) 

25.9 

25.9 

25.9 

4,  (s-1) 

100.53 

io12-91 

io12-21 

42  (s-1) 

io1778 

1018-14 

IO1790 

Cl 

0.36 

0.31 

0.26 

C2 

0.35 

0.36 

0.47 

where  the  weight  factor  c,  is  equal  to  the  amount  of  volatiles  formed 
from  the  ith  pseudocomponent,  and  M  is  the  number  of 
pseudocomponents. 

It  is  noted  that  a  pseudocomponent  represents  the  totality  of 
those  decomposing  species  that  can  be  described  by  the  same 


kinetic  parameters.  Those  pseudocomponents  are  usually  linked 
the  real  biomass  components.  In  particular,  for  the  three-parallel- 
DAEM-reaction  model  [25,55],  the  first  and  second  pseudocompo¬ 
nents  represent  the  fractions  of  hemicellulose  and  cellulose  that 
react  at  low  temperature,  and  the  third  pseudocomponent  repre¬ 
sents  the  sum  of  the  remaining  amounts  of  hemicellulose  and 
cellulose  plus  the  fractions  of  lignin. 

By  following  the  above  stepwise  discussion,  the  research  and 
development  efforts  in  the  literature  on  the  DAEM  for  lignocellu¬ 
losic  biomass  pyrolysis  are  summarized  in  Table  4.  Most  kinetic 
studies  of  lignocellulosic  biomass  pyrolysis  have  been  made  using 
thermogravimetric  analysis  (TGA)  at  relatively  slow  heating  rates. 
In  the  tabulated  summary,  the  elemental  composition  of  lignocel¬ 
lulosic  biomass  sample,  the  experimental  methods  and  conditions 
are  highlighted.  The  resultant  DAEM  kinetic  parameters  together 
with  the  parameter  estimation  methods  are  also  provided.  The 
efforts  on  DAEM  for  the  pyrolysis  of  lignocellulosic  biomass 
indicate  that  the  model  consisting  of  three  parallel  DAEM 


J.  Cai  et  al.  /  Renewable  and  Sustainable  Energy  Reviews  36  (2014)  236-246 


245 


Table  4.2 

Kinetic  parameters  of  lignocellulosic  biomass  pyrolysis  in  the  literature  [25], 


Corn  stalk 

Sorghum  straw 

Rice  husk 

Wheat  straw 

Eoi  (kj  moL1) 

176 

176 

176 

176 

Eo2  (kj  moD1) 

185 

185 

185 

185 

Eos  (kj  mol-1) 

189 

189 

189 

189 

ax  (kj  mol-1) 

7.1 

7.1 

7.1 

7.1 

0-2  (kj  mol-1) 

1.7 

1.7 

1.7 

1.7 

o-3  (kj  mol-1) 

32.7 

32.7 

32.7 

32.7 

Ai  (s-1) 

10™ 

101413 

101413 

101413 

A2  (s-1) 

101371 

io13-71 

1013-71 

101371 

43(s— ’) 

10139° 

O 

o 

io13-90 

10139° 

Cl 

0.14 

0.10 

0.08 

0.16 

C2 

0.25 

0.23 

0.36 

0.38 

C3 

0.31 

0.34 

0.33 

0.19 

Table  4.3 

Kinetic  parameters  of  lignocellulosic  biomass  pyrolysis  in  the  literature  [55], 


Corn  stover 

Cotton  stalk 

Palm  oil  husk 

Pine  wood 

Red  oak 

Sugarcane  bagasse 

Switchgrass 

Wheat  straw 

£01  (kj  mol-1) 

179.602 

178.188 

169.710 

186.701 

183.109 

184.750 

186.776 

175.506 

£02  (1<J  moU1) 

207.381 

205.287 

199.966 

214.364 

209.58 

212.483 

212.057 

204.244 

Eos  (1<J  moU1) 

239.344 

239.463 

236.110 

271.758 

242.145 

234.759 

260.952 

240.610 

o i  (kj  mol-1) 

5.888 

5.591 

5.147 

8.769 

8.342 

5.419 

4.586 

5.375 

a 2  (kj  mol1) 

1.174 

1.789 

1.313 

1.126 

0.706 

1.339 

1.361 

1.130 

a 3  (kj  mol-') 

31.414 

41.767 

39.998 

29.250 

26.583 

36.493 

39.293 

38.422 

A,  (s-1) 

|Q  13.007 

1013132 

|q12.699 

|q13.006 

|q13.007 

|q13.007 

|q13.007 

|Q  12.700 

A2  (s-1) 

|Q  13.949 

|q13.800 

|q13.618 

|q13.803 

|q13.618 

|q13.805 

|q13.804 

|q13.800 

A3  (s-1) 

|q16.003 

io15-738 

io15-953 

io16-399 

|q15.006 

io15-769 

|q16.536 

|q15.812 

Cl 

0.16337 

0.12999 

0.10116 

0.32147 

0.27282 

0.22009 

0.23651 

0.15916 

c2 

0.25386 

0.35824 

0.32121 

0.38557 

0.37320 

0.35035 

0.33376 

0.28027 

C3 

0.16481 

0.20008 

0.19241 

0.08412 

0.12204 

0.12289 

0.20159 

0.23410 

reactions  is  the  most  accurate  and  up-to-date  approach  for 
modeling  the  pyrolysis  kinetics  of  lignocellulosic  biomass. 

7.  Summaries 

The  pyrolysis  of  lignocellulosic  biomass  involves  many  reac¬ 
tions,  which  can't  be  well  described  by  a  simple  single  kinetic 
model.  In  the  literature,  the  DAEM  was  usually  used  to  describe  its 
pyrolysis  kinetics.  This  review  presents  a  state-of-the-art  review  of 
the  DAEM  for  lignocellulosic  biomass  pyrolysis.  This  review  starts 
with  the  derivation  of  the  DAEM  and  a  brief  overview  of  the 
activation  energy  distribution  and  frequency  factor  in  the  DAEM. 
Then  the  parameter  estimation  methods  for  the  determination  of 
the  DAEM  kinetic  parameters  are  critically  studied.  The  final  part  is 
the  overview  of  the  application  of  the  DAEM  to  lignocellulosic 
biomass  pyrolysis.  A  handy,  up-to-date  overview  of  the  existing 
research  and  developments  efforts  and  major  findings  in  the 
DAEM  for  lignocellulosic  biomass  pyrolysis  is  provided.  It  is  hoped 
that  this  review  is  critical  for  helping  understand  the  application 
of  the  DAEM  in  the  thermal  decomposition  of  other  degradable 
materials. 
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